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Phase transitions involving spontaneous time-reversal symmetry breaking are studied on the hon¬ 
eycomb lattice at finite hole-doping with next-nearest-neighbor repulsion. We derive an exact ex¬ 
pression for the mean-field equation of state in closed form, valid at temperatures much less than 
the Fermi energy. Contrary to standard expectations, we find that thermally induced intraband 
particle-hole excitations can create and stabilize a uniform metallic phase with broken time-reversal 
symmetry as the temperature is raised in a region where the groundstate is a trivial metal. 
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Introduction .— A popular proposal to break time- 
reversal (T) symmetry in an electronic system without 
invoking the spin degrees of freedom and without break¬ 
ing translational symmetry is to imagine microscopic 
electronic current loops arranged so that their net mo¬ 
ment vanishes in each unit cell [1]. In multiband sys¬ 
tems, T symmetry breaking is intimately related to the 
condensation of interband particle-hole pairs when the 
electron-electron scattering is strong enough [2]. The re¬ 
duced symmetry of the vertex function Dk that couples 
the interband bilinear operators a^bk determines both 
the symmetry breaking and the symmetry of the order 
parameter ip ~ In particular, T symmetry 

will be broken if T'DkT~^ = TA*_k ^ [3-5]. 

Metallic systems with non-trivial topological proper¬ 
ties, generically referred to as topological Fermi liquids 
(TFLs) [6], have come to play an important role in the 
study of marginal and other non-Fermi liquid phenom¬ 
ena [7]. Various instances [8-16] of TFLs and topologi¬ 
cal Mott insulators (TMIs) at commensurate fillings have 
been realized on the honeycomb lattice using a combina¬ 
tion of mean-field theory and numerics. 

Although the breaking of T symmetry opens a gap at 
the Dirac points [9], no gap opens at the Fermi surface 
in mean-field theory [5]. Thus a continuum of gapless 
intraband particle-hole excitations exists in the ordered 
phase in one-to-one correspondence with the free Fermi 
gas. They provide the leading temperature corrections 
to the free energy A) at the mean-field level. 

Here, /r is the chemical potential, A is the external field 
that couples to the order parameter, and T is temper¬ 
ature. The order parameter ip is given by the relation 
dUlyi-p/dX = —ip] the resulting integral equation is in 
general not restricted to the vicinity of the Fermi sur¬ 
face but extends over a large portion of the Fermi sea 
up to a cut-off scale A ^ 1/a, where a is some micro¬ 
scopic scale determined by the range of the interaction. 
This results in an increased sensitivity to the number of 
particles N within the cut-off energy. Thus the equation 


for the order parameter has to be solved simultaneously 
with the number equation = —N. The solu¬ 

tions, if derivable, can be used to construct the canonical 
free energy Fmf(T N, (p) = Dmf(T F, A) -b -b pN. 

In this paper, we obtain an analytical solution for 
the mean-field equation of state dFMF/dp = ^ and 
dFur/dN = p ior N interacting spinless particles with 
next-nearest-neighbor repulsion on the honeycomb lat¬ 
tice away from half-filling. The solutions are valid in the 
limit T ep, where is the Fermi energy; physically, 
in this limit only the intraband particle-hole excitations 
around the Fermi surface provide a significant contribu¬ 
tion to the partition function at finite T, while the pair¬ 
breaking interband excitations can be neglected. 

Analysis of these equations in the limit A —0 allows 
us to identify and study both the quantum critical point 
(QCP) and the low temperature phase diagram. We find 
that the thermal smearing of the Fermi surface has the 
curious effect of favoring the condensation of interband 
particle-hole pairs resulting in a finite temperature phase 
transition into an ordered state with broken T symmetry. 
This leads us to the novel conclusion that a stable topo¬ 
logical state can be created by raising the temperature. 

Model .— We consider the extended Hubbard model of 
spinless electrons on a honeycomb lattice with nearest- 
neighbor [nn) hopping and next-nearest-neighbor (nnn) 
repulsion. It was concluded in Ref. [15], after extensive 
mean-field studies of this model at T = 0, that a finite 
nnn repulsion is necessary to stabilize a T symmetry 
broken metallic phase in the vicinity of half-filling, which 
is consistent with Ref. [9] where the half-filled case was 
first studied. They compete with other charge-modulated 
states that occupy most of the parameter space when nn 
repulsion is included. Here, we isolate and study the 
stability of the T broken phase at finite T, keeping only 
the nnn repulsion while neglecting the nn repulsion. 

For the Cqv point symmetry group of the underlying 
honeycomb lattice, only the one-dimensional irrep (called 


2 


Bi [17]) with a reduced Cq symmetry 


Bk 


1 

7 ! 


[sin(fc • ti) + sin(fc • t^) + sin(fe • t^)] 


( 1 ) 


is real and odd so that ^ Vk [5]. (The vectors ti, 
with I = 1,2, and 3, are the basis vectors of the hexagonal 
Bravais lattice.) The effective Hamiltonian that isolates 
this irrep in the particle-hole sector reduces to [5] 

^ = ^'^k[^kT:3- ( 2 ) 

k q 


( ^k,k' '^^k — k''^kj^fk^k'- TllG 

shorthand notations, k = {k, ikn) and k — k' = q = 
(Q,i<Zm), where fc„ (q-m) are the odd (even) Matsubara 
frequencies, are used; the Fourier transform is defined as 
'/'(a;) = The operator fk,k' 

represents all the terms within the square brackets in (3); 
note fk,k = 'Ti- The momentum = {k + k')/2 in Vk^. 

Assuming a uniform solution (pq = /3Sm ,hr S, 
followed by a few standard manipulations involving the 
identity Tr In A = In det A and summing over the fre¬ 
quencies [21, 22], the action after scaling S = {I3L"^)C 
can be written as £ = ip‘^l{2V) + £+ + £_, where [14] 


The basis tpl, = [b],, a7/72 diagonalizes the kinetic term 
after a unitary rotation ak = e~^^’°Ak+e^^>°Bk and bk = 
g-i^hAk — e^^'^Bk, where Ak and Bk are the momen¬ 
tum space representations of the on-site annihilation op¬ 
erators on the two sub-lattices of the honeycomb lattice. 
The Ti’s are Pauli matrices that act on the band indices. 
The angles are obtained from exp(zfe • di), 

where the vectors di connect the nn atoms on the two 
sub-lattices. The band energies €k = '/3~+~2Ck, where 
Ck = cos(fc • ti), have vanishing gaps at two inequiv¬ 
alent Dirac points, ±ku. All energies are measured from 
the Dirac points in units of the hopping integral t. 

The bilinear operator in the interaction term reads [5] 


$(q) 




k+i 


r(efc+3-efc_,)l 




V'fc-2 (3) 


Although our choice of the ip basis complicates $ by in¬ 
troducing additional phases, defining the order parame¬ 
ter a.s (fi ^ {^{q — 0)) makes explicit the mapping of K 
in (2) to the transverse-field Ising model: The system un¬ 
dergoes a transition at a critical Vc from a “paramagnet” 
with all spins pointing along T3 = — 1 (fully occupied 
lower band when (p = 0) to a correlated “Ising ferro- 
magnet” with the spins pointing along Ti when tp ^ 0. 
It is well known that the collective excitations of the 
transverse-field Ising model are gapped [18, 19]. How¬ 
ever, in a metallic system the low-energy fermionic exci¬ 
tations qualitatively changes the nature of the transition 
[5, 14, 20]. This is the key point that we analyze here. 

Fixed pi. — The action corresponding to K\ = K — 
X^{q = 0), where the Hamiltonian K is defined in (2) 
and A is a constant source term that couples to the field 
operator $(<7 = 0 ) = J2k j is derived the stan¬ 

dard way [ 21 ] by first decoupling the quartic fermion in¬ 
teraction terms using Hubbard-Stratonovich (HS) fields. 
Since $ is Hermitian, the HS fields, (pq, are real, i.e., 
p* = p-q', note that, when q — 0, the phase factors in 
(3) drop out leading to considerable simplifications. 

The HS transformation allows the fermions to be for¬ 
mally integrated out. The resulting partition function is 
written as Z = f D[p] exp(—S'), where the bosonic action 
[5] S = (TVW)Eq </>-,'^9 - TrlnG-1 with = 


£± = ^In (1 + e (4) 

k 

The renormalized bands equal = FEx k — Mi where 
E\k = (cfc + p\ = p + \ is the shifted order 

parameter and the bare energies are defined in ( 2 ). 
From here on, we confine ourselves to the linear regime 
around the Dirac points. As is well known, the bare 
spectrum is linear Cp = a\p\ to first order in p — k — kjj 
with a = 3/2. Furthermore, we assume 
note that coincidentally the constant = 3/2 = a, 

as a result the band energies assume the simpler form 

+ + (5) 


The linear approximation requires a cut-off momentum 
PA, which we determine by fixing the number of electrons 
in the cones at T = 0. It is easily seen that the mean- 
field Hamiltonian [HMF,f>-k] = 0 conserves the particle 
number in each k [5], hence the number of electrons Ng 
is always equal to the number of states between the non¬ 
interacting cut-offs pf and pA independently of p, i.e.. 



A - 9 

V L 2 - 


^(Pa-Pf) ( 6 ) 


The factor 2 accounts for the two Dirac cones. For 
definiteness, we consider hole doping; the hole density 
rik =n — ne, where n = e^/(27ra^) and = P^p/(2'Ko?'\ 
and the bare energies ca = apA and ep = app- It is 
convenient to write the renormalized cut-offs in terms of 
the rescaled densities n = (Tra'^)n and he = {Tra^)ne as 
Aa = a/ca T = x/2h + A| and Ex^f = = 

x/2.hh + A^. It follows that 2|Aa| = 2a\px\ is the gap at 
the Dirac points. Eq. ( 6 ) then reduces to 

he = h-hh = ]^{Al-EI p) (7) 

which is independent of Aa as discussed. 

The terms C± defined in Eq. (4) are valid for all tem¬ 
peratures and fillings. Here, we restrict ourselves to the 
temperature range T/Ex^f 1 , which corresponds to 
suppressing all terms that fall off exponentially as 
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and higher. (This approximation breaks down at half¬ 
filling, where /r = 0 from particle-hole symmetry. For a 
detailed analysis of the half-filled case, see [23].) 

To this end, we first evaluate the momentum sums 
in (4), which can be done analytically, and then ex¬ 
panded using the asymptotic properties of the polyloga¬ 
rithm functions [24]. The result of the expansion when 
combined with the l{2V) term in C gives (see appendix 
[Al] for the derivation of Eq. ( 8 )) 

- A3 Ail/r] 1 m|3 

+ ^ 6 ^ 

Here, C = ( 7 ra^)£ (similar to the rescaling of n). Since 
hole doping is assumed, /r = — J/r] 0. The contribu¬ 
tion is the familiar Sommerfeld correction that originates 
from the particle-hole excitations at the Fermi surface; it 
plays a key role in what follows. 

The mean-field free energy r2(/i. A) is defined as 
H(/i, A) = £(A,/r, A), where A satisfies the saddle-point 
equation 9£/9Aj^ = 0 [25]; differentiating ( 8 ), we get 

TT /\ — — 

— =A,(A,-H) (9) 


(The bars denote that the quantities are evaluated at the 
saddle-point.) The same equation is obtained by setting 
the total derivative dC/d\ = —(7ra^/L^)(4(q = 0)) = 
—a^irA/V) at the saddle-point, which is consistent with 
assigning A as the order parameter in the limit A —^ 0. 

Fixed Tie. — Since H as defined is the free energy, the 
particle density corresponds to he = differen¬ 

tiating (8) again, we get he = (A| — /r^)/2 — 7r^T^/6. 
Comparing he with its T = 0 value in (7), we obtain 

^x{Tf = EIp - -T^ = 2hh + Al-— (10) 

Solving Eqs. (9) and (10) for A and fi in terms of A 
and he and substituting the solutions in 

- TT A 

F{A,he) = id{n,X) + aX—+ nhe (11) 
gives the free energy F{A,he) (per unit volume). 


Fortunately, it is possible to solve (9) and (10) in closed form. For the equation of state, we get (See appendix [A2] 
for the derivation of Eq. (12)) 


aA = 



1 / *? 

/ {2h + 2hh - tt^TVO) -f 2^2h{2hh - tt‘^T'^1^) + (7rA/y)2\ V 

y (2h — 2hfi + — 4(7rA/V)^ J tt 


( 12 ) 


Substituting Eq. (12) back into (10) gives the solution for j^j. Note that the cut-offs are absorbed in h and h^. 


r 


At r = 0 , the QCP, Vc, is obtained by setting A and 
A = 0 in (12). This reproduces the result [5] 


Ve = 


— \/2hh CA — ep 


(13) 


Our key observation is that when T > 0, a non-trivial 
solution emerges in the region V < 14 at a finite temper¬ 
ature Te- To see this, we set A = 0 in (12) and get 


V{Te) 


TT 

A/^-7r2T2/3 




(14) 


Eq. (14) can be inverted to give the finite temperature 
phase boundary. For V close to 14, we get 


a, r ^ n- (15) 

This is an intriguing result: contrary to standard expec¬ 
tations, a topologically non-trivial metallic phase with 
broken F symmetry, which is unstable at T = 0 for 
V < 14 , is stabilized as the temperature is raised above 
Te- The prefactor is a reminder that the term 


originates from the intraband particle-hole excitations; 
hence, this thermally induced finite temperature T sym¬ 
metry broken phase only exists at finite doping. 

It is easily verified that Eq. (14) can be obtained di¬ 
rectly from (9) and (10), instead of using the solution in 
(12). Neither way, however, guarantees that the A 7 ^ 0 
phase is stable and is lower in free energy than the trivial 
phase with A = 0 . For this, the difference in the free en¬ 
ergies SF{A) = F{A,he) — F{0,he) has to be analyzed. 

Although it is straightforward to calculate SF given 
the equation of state ( 12 ), the algebra is tedious and the 
result in its full generality is not very informative. The 
fact that the state is stable can easily be seen by plotting 
it numerically as shown in Fig. 1. A clear minimum is 
observed when T > Te in the region where E < 14- 

The curvature at the minimum can be obtained ana¬ 
lytically by differentiating the equation of state in ( 12 ), 
{aTT/V)X = dT/dA\n^, holding he fixed: 


d'^F 


aA 2 


OTT dX 


aA^ 


> 0 as T ^ Te (16) 


The derivative is evaluated at E = E (T[4), slightly above 
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FIG. 1. The difference of the canonical free energy (5J^(A) 
defined in (11) is evaluated numerically keeping the density 
he fixed and plotted for V < 14 as a function of A for T = Tc 
and T > Tc- A clear minimum is observed as the temperature 
is raised above Tc (see Eq. (15)), signifying the formation 
of a topological liquid phase at finite T. The evolution of 
the order parameter A(r) (the location of the minimum) is 
shown in the inset by the solid line; the dashed line is the 
approximate expression for A derived in (19). The filling is 
fixed at = •Jn^jn = 0.1 and 14 = 1.3. 


the critical temperature where A is small but finite. 

Approximate analytic expressions for the variation of 
the order parameter A as a function of 14 and T when 
1 are easily obtained by expanding ( 12 ) 
near V),. Setting A = 0 and expanding (12), we get 


CACF CF V{TY \ V{TY ) 


(17) 


Here epiT) = \/£% — and the function 14(T) has 

the same form as V{Tc) in (14) with Tc replaced by T. 
Near Tc, V(T) can be expanded to linear order in (Tc/€f) 
in powers ot ST = T - Tc as 


VjTc) - V{T) _ F 
V{Tc)^ ^ 3 


{STY 


^)ST-. 
ep J 2eF 


(18) 


From (14), we see that Tc —?► 0 as V{Tc) —>■ 14) hence 
Tc 0 only in the weak-coupling regime 14 < I 4 . The 
(ST)^ term is therefore necessary in the strong-coupling 
regime V > Vc as the linear term vanishes there. 

The behavior of A as 14 crosses from V~ to is 
analyzed below using Eqs. (17) and (18): 

{i) For 14 < Vc, we fix the potential at 14 = V{Tc) in 
(17) and expand V (T) = V (Tc+ST) using (18) to leading 
order in ST/Tc and get 


A(T)2 = 


27r^T? fT-T 


as 14 —>• 14 


(19) 


A plot of this function indicated by dashed lines in the 
inset to Fig. 1 agrees accurately with the numerically 
extracted minimum (solid line in the inset) as T —> T^. 


{ii) For 14 > 14 , we set Tc = 0 and 14 = Vc + SV in 
(17). From (18) weget E(T)/14 « ^-{.T^Vcep/&){T/epY, 
which when substituted in (17) gives 

A(T)2 = A(0)2 + ^ as E ^ 14 + (20) 


where A( 0 )^ = {2Ti){eF/Vc){SV/Vc) gives the standard 
mean-field exponent of 1/2 at T = 0 [5]. 

{in) Finally, at 14 = I 4 , since A(0) = 0 at the QCP, 
we get from (20) that A(T) = FT/\fi\ the linear T de¬ 
pendence is distinct from that on either side of the QCP. 

To conclude, our main observation from Eqs. (19) and 
( 20 ) is that contrary to standard expectations, the in¬ 
traband excitations at finite T has a stabilizing effect 
that tend to increase A as T increases. The increase 
will eventually be limited by the interband excitations 
when T > ep, the study of which is beyond the scope 
of this paper. More precisely, we note from Eq. (15) 
that Tc/ep ^ \/(14 — V)/eF- Hence, provided <514 ^ ep, 
pair-breaking effects are parametrically far and our re¬ 
sults are internally consistent with our approximations. 

Finally, we want to stress the importance of studying 
the canonical free energy T{ne) (since the volume is held 
constant throughout, this translates to fixing the den¬ 
sity), rather than directly minimizing the action T(/i) 
holding /j, constant. Although demanding dC/dA — 0 
recovers the saddle-point condition in (9), the conden¬ 
sate is unstable. To see this, we differentiate ( 8 ) twice 
and obtain d'^C/dA'^ = —A^/A < 0. This instability 
was pointed out earlier by us in Ref. [5] by an explicit 
diagrammatic calculation. In Eq. (16) we resolve this in¬ 
stability by showing that the mean-field condensate ob¬ 
tained keeping N fixed, rather than /r, is indeed stable. 

Notice that the order parameter preserves inversion (T) 
symmetry but breaks chiral (C) symmetry in such a way 
that the combined C'T symmetry is left invariant, as a re¬ 
sult it can exhibit both anomalous Hall and Kerr effects 
[4, 6 , 8 ]. Furthermore, as noted in Ref. [4], because of 
the chirality, the state does not couple directly to non¬ 
magnetic impurities. Hence when spin effects can be ig¬ 
nored, the phase is expected to be stable. It remains to be 
seen, however, if the finite temperature phase described 
in this work survives in the presence of nearest-neighbor 
repulsion, which we have not included in this study. 

Temperature induced topological transition has been 
reported in the context of thermally induced band in¬ 
version due to electron-phonon interactions [26]. In our 
case, the topological state is induced by strong correla¬ 
tions and thus provides a qualitatively different paradigm 
involving only electronic degrees of freedom. 

AP acknowledges E. Ponton (IFT), J. Birman and M. 
C. N. Fiolhais for helpful discussions, and we thank P. 
Ghaemi for drawing our attention to Ref. [26]. 
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APPENDIX 


lipses indicate higher powers of e 


[Al] DERIVATION OF THE LOW 
TEMPERATURE LIMIT OF £ 


In this section, we detail the steps taken to arrive at 
the low temperature limit for the action C derived in Eq. 
( 8 ). We start from the definition 

£ = (/?V( 2 U) + £++ £_ 


where C± are defined in Eq. (4). Eollowing the discus¬ 
sions leading up to Eq. ( 8 ), it can be shown that the 
sum over k in C± can be expressed as an integral with 
the appropriate cut-offs as 

£± = ^ ln(l + 

where /i = — |/r|. Substituting z± = (ay^ + £^±H)/t 
brings the integrals into the more familiar form: 

^ (Aa±|m|)/T 

C± = —^ [ dz(rz± =F lMl)ln(l 

TTa^ J 

{\A^\±M)/T 


The above integrals can be represented in terms of the 
polylogarithm functions Li 2 and Li 3 . (See Ref. [24] for 
the definitions and properties of the polylog functions.) 
Rescaling £ = {Tra‘^)£, we obtain 


£+ = 


Ai |A,|3 , Al\^,\ All 


+T^ 

_ rp3 


AT* / ^A + I^l \ I . It- . / I^aI + I^I \ 

AaLi 2 (^-e T j - |Aa|Li 2 (^-e ^ j 


Li .3 (-e T j - L 13 (-e T ] 


£- =T^ 


AaLi 2 (^-e T j _ |Aa|Li 2 (^-e j 


_ 'J^O 


-e T 

AaI-I>^I 


T. / \ / |AaI-I>x| a 

L 13 (-e T j - L 13 (-e T ] 


lim Li 2 (—e 

Re(2)—>-oo 

Li 2 (-e^) 

Li3(-e"^) 

Li 3 (-e^) 




H- 


—e 




Dropping terms that fall off as and higher, we 

get £+ = 0. Note that |Aa| < Ea.f = a/A^ -I- away 
from half-filling and A a ^ e\. Gathering the polynomial 
terms from £_ and substituting into £ gives Eq. ( 8 ). 


[A2] DERIVATION OF THE EQUATION OF 
STATE 


In this section, the algebraic steps leading up to the 
solution for the equation of state derived in ( 12 ) are de¬ 
tailed. This requires solving for A and /x using the saddle- 
point equation (9) and the equation relating the chemical 
potential and density ( 10 ). 

The relevant equations to be solved are listed below 

■jr /\ _ _ 

Eq. (9) : — = Aa (Aa - l/xj) 

7r2y2 

Eq. (10) : = 2nh -f - — 

where Aa = V2n + A^. Our aim is to eliminate the 
square-roots from Aa and /x and obtain a polynomial 
equation that can be easily solved for Aa. To this end, 
we multiply both sides of the first equation by (Aa + |ax|), 
and use the second equation to eliminate /x^. This sim¬ 
plifies the equations to 


Aa + ImI 
Aa - ImI 



2fih + 7r2rV3)AA 


To obtain the low temperature expansion, we use the 
asymptotic properties of the polylog functions. The re¬ 
quired formulas to order e~^ are listed below: (the el- 


It is now straightforward to solve for Aa by adding 
the two equations and squaring the result to eliminate 
the square-root. After some minimal algebra, we get 


/ttA^ / (2fi -b 2fih - tt^T'^/ 3) + 2^2fi{2fih - Tr^T^/i) + (ttA/V)^ 

V“i^ J (2fi - 2fih + 7r2T2/3)2 - 4(7rA/U)2 


Care must be taken in choosing the correct sign of the equation for A|. It is chosen to reproduce the solution 
square-root in the numerator when solving the quadratic for V), in Eq. (13), which can be independently obtained 
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from the saddle-point equation in the limit A = A = 0. 

Finally, we take the positive root of A| to obtain the 
equation of state given in (12). This is followed by sub¬ 
stituting the solution in the equation for to find |p,|. 
The solutions for A and |/i| are used to obtain the exact 
low-temperature canonical free energy JF^ defined in Eq. 
(11). The result for bT is plotted in Fig. 1. 
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